QUANTUM HAMILTON - JACOBI STUDY OF WAVE FUNCTIONS 
AND ENERGY SPECTRUM OF SOLVABLE AND QUASI - 
EXACTLY SOLVABLE MODELS 




SCHOOL OF PHYSICS 
UNIVERSITY OF HYDERABAD 
HYDERABAD - 500046, INDIA 
DECEMBER 2003 



1 



Declaration 



I, K. G. Geojo, hereby declare that the work reported in this dissertation titled, 
Quantum Hamilton - Jacobi study of wave functions and energy spec- 
trum of solvable and quasi - exactly solvable models, is entirely original 
and has been carried out by me, under the supervision of Prof. A. K. Kapoor, 
Department of Physics, School of Physics, University of Hyderabad. 

To the best of my knowledge, no part of this dissertation was submitted for any 
degree of any other institute or university. 



Place: Hyderabad 

Date: K. G. Geojo 



Certificate 



2 



This is to certify that the report entitled Quantum Hamilton - Jacobi study 
of wave functions and energy spectrum of solvable and quasi - exactly 
solvable models, being submitted by K. G. Geojo, in partial fulfillment of the 
requirements for the award of Doctor of Philosophy in Physics by University 
of Hyderabad, is a bonafide work carried out at the University of Hyderabad 
under my supervision. The matter embodied in this report has not been submitted 
to any other institute or university for the award of any degree. 



Dean, 

School of Physics, 
University of Hyderabad, 
Hyderabad - 500046. 



Prof. A. K. Kapoor, 

School of Physics, 
University of Hyderabad, 
Hyderabad - 500046. 



Contents 



1 INTRODUCTION 



2 QUANTUM HAMTLTON-.T ACORT FORMALISM 




2.1 Classical Hamilton- Jacobi Theory 


2.2 Quantum Hamilton- Jacobi Eauation 




2.3 Exact Quantization 






2.3.1 Boundary condition 


2.3.2 Exact Quantization Condition 






2.4 Connection with Schrodineer Eauatio 


n 


2.5 Singularities of QMM 


2.6 Energy Spectrum of Morse Oscillator 





3 CACULATION OF WAVE-FUNCTION FOR ES MODELS 



3.1 Harmonic Oscillator 



3.2 Morse Oscillator! ■ ■ ■ ■ 

3.3 Poschl- Teller Potential 



3.4 Eckart Potentia 



3.5 Hydrogen Atom 



4 CONDITIONS FOR QU AST-EX ACT SOLVABILITY 




4.1 Introduction to OES 





CONTENTS 



4 



4.5 Circular Potential! 45 



4.2 A Representation of QES Quantization Ruld 40 



4.3 Sextic Oscillator 41 



4.4 Sextic Oscillator with a Centrifugal Barrier! 43 



4.6 Hyperbolic Potential! 4 



4.7 V(x) = A sinh 2 \ fi)x + B sinh ^fvx + C tanh -t/t/x secfu/z/x + ,Dsech 2 -t/z/x 

. _ , . . O , ~ . , ^ , , , _ , O , 



52 



4.8 V(x) = A cosh 2 a/zae + B cosh 4- C coth ^/z/xcsch^/^'a: 4- -Dcsch^/z/x I 55 



4.9 V(x) = Ae 2 ^ x 4- Be^El + Ce"^ + £>e~ 2 ^ 



58 



4.10 Quartic Oscillator! 61 



4.11 Summary and Observation; 



62 



5 


CALCULATION OF WAVE-FUNCTIONS FOR OES MODELS 




5.1 Sextic Oscillator 






5.2 Hyperbolic Pote 


ntia 






5.3 Concluding- Remarks 





69 
74 



6 CONCLUSIONS AND OUTLOOK 



76 



Chapter 1 
INTRODUCTION 



In this thesis we present an alternative approach to the study of exactly solvable 
and quasi-exactly solvable (QES) problems in quantum mechanics. This approach, 
known as quantum Hamilton- Jacobi (QHJ) approach, [1] has been found to be an 
elegant and simple method to determine the energy spectrum of exactly solvable 
models in quantum mechanics. The advantage of this method is that it is possible to 
determine the energy eigen-values without having to solve for the eigen-functions. In 
this formalism, a quantum analog of classical action angle variables [2] is introduced. 
An exact quantization condition is formulated as a contour integral, representing the 
quantum action variable, in the complex plane. This exact quantization condition 
has been utilized for determining the energy eigen-values for one dimensional and 
separable systems. The quantization condition represents well known results on the 
number of nodes of the wave-function, translated in terms of logarithmic derivative, 
also called quantum momentum function (QMF). The equation satisfied by the QMF 
is a non-linear differential equation, called quantum Hamilton- Jacobi equation leads 
to two solutions. A boundary condition — in the limit h — > QMF tends to the 
classical momentum — is used to determine physically acceptable solutions for the 
QMF. The application of QHJ to eigen-values has been explored in great detail by 
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Bhalla et al [3,4]. 

In chapter 2 we review QHJ method and, by means of an example, we show how 
eigen-values are calculated without the need to obtain the full wave function. Briefly, 
this is possible because for the implementation of exact quantization condition one 
needs the knowledge of the singularities of QMF and the residues. The residues are 
easily computed by substituting only a few terms of the Laurent's expansion in the 
QHJ equation. 

In chapter 3 we show how to calculate bound state wave functions in the QHJ 
formalism. For this purpose, again, one only requires knowledge of singularities of 
QMF and the corresponding residues. The technique to calculate residues is already 
available from earlier works and these are used for obtaining the bound state wave 
functions. In this process we clarify certain assumptions which are needed, and 
are found to be correct, for all the exactly solvable models which we have studied. 
As a by-product of the study of the bound state wave functions we have another 
way of obtaining the energy eigen-values. In this chapter we present details of our 
calculation for harmonic oscillator, hydrogen atom, Poschl Teller, Morse and Eckart 
Potentials, while details of some other potentials can be found in [5]. 

In chapter 4 we take up a study of the QES potential models [6] in one-dimension. 
These models have been extensively studied using Lie algebras. The QES mod- 
els have the property that a part of the energy spectrum and corresponding wave- 
functions can be computed exactly if the potential parameters satisfy a constraint 
known as the condition for quasi exact solvability. We study several QES models 
within the frame work of QHJ method. In each case we show that the condition for 
quasi exactly solvability follows from a very simple assumption about the behavior 
of QMF at infinity. Our assumption is equivalent to assuming that, after a suitable 
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transformation, the QMF reduces to a rational function of the independent variable. 
For all the known QES models, the quasi-exact solvability condition can be derived 
in this fashion [7]. 

In chapter 5 we study the wave functions of quasi-exactly solvable models and 
present details of our calculation for sextic oscillator and hyperbolic potential. We 
find that obtaining eigen-values and eigen-functions does not require any new tech- 
nique other than those given in chapter 3 for the exactly solvable models. However, 
this study reveals an interesting property of QMF for the bound states of quasi- 
exactly solvable models. This result concerns the zeros of the bound state wave 
functions in the complex plane. In the case of exactly solvable models, the moving 
poles of QMF appear only on the real line and all such poles correspond to the nodes 
of the wave- function. The number of such poles increases with energy in accordance 
with well known theorems on the number nodes of the wave-function. In the case of 
QES models all the bound state wave-functions, which are computable algebraically 
and also by our method, have complex zeros in addition to the real zeros correspond- 
ing to the nodes. In fact we find that for a given QES potential all such 
wave-functions have the same number of zeros if we count all real and 
complex zeros. 

In the last chapter, we give a summary of our work as well as some directions for 
further investigations within the QHJ formalism. 



Chapter 2 



QUANTUM 
HAMILTON- J ACOBI 
FORMALISM 

In this chapter, we summarize the main results of the Hamilton- Jacobi theory in 
classical mechanics and the QHJ formalism to be used in this thesis. In section 
3 a quantization condition is given, which is exact for one dimensional system and 
separable systems in higher dimensions. In section 4 the connection of QHJ formalism 
and Schrodinger quantum mechanics is spelled out and in the last section of this 
chapter an example of computation of eigen-values within the QHJ formalism, and 
without solving for wave-functions, is given. 

2.1 Classical Hamilton- Jacobi Theory 

The phase-space formalism of classical mechanics gives us freedom to introduce a 
pair of cannonical variables (Qk, Pk) which are functions of a given starting set of 
variables (qk,Pk)- The Hamiltonian form of equations of motion is preserved if the 
transformation is cannonical in the sense of preserving Poisson brackets. This free- 
dom is utilized in the Hamilton- Jacobi theory to give a formal solution of classical 
mechanical problems by making a transformation, so that the Hamiltonian becomes 

8 
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constant. 

In the Hamilton- Jacobi theory we look for a function W(q i: Pi), which generates 
the desired cannonical transformation making the Hamiltonian a constant. The 
transformation equation relating the old and new cannonical variables are 

dW „ dW 

p -=%-* = m (21) 

and the requirement, that the Hamiltonian in terms of new variables be a constant 
cti, gives a partial differential equation for W(qi,Pi): 

H {*•!&) =ai - (2 ' 2) 

This equation is the Hamilton-Jacobi equation. The function W(qi, Pi) is known 
as the Hamilton's characteristic function. It is well known that a solution to the 
Hamilton-Jacobi equation is equivalent to full solution of Euler Lagrange equations 
of motion [2]. 

In general W(qi, Pi), which is a function of q; L and Pi, can be taken to be a function 
of (?i's and constants of motion ol\, a% • • • , by identifying the new momenta Pi with 
the constants of motion a iy where ot\ is the total energy of the Hamiltonian. 

For the purpose of finding the frequencies without solving the equation of motion 
completely, action variable Jj is introduced by 

Ji={ dW ^ a ^'--^ dqi (2.3) 
J dqi 

where the integral is over a periodic orbit. The action variable Jj are functions of 
the constants of motion a, and one can eliminate a^s in favor of the action variables 
Jj. In particular the Hamiltonian, H — a±, can now be expressed in terms of action 
variable J« 

H = H(J 1 ,J 2 ,---,J n ) (2.4) 
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The generalized phase-space variable conjugate to Jj are known as the angle 
variable u>i and are given by the transformation equations 

9H(Ji, J 2 , ■ • ■ , J n ) , . , . 

w<= -qJ = v i (J 1 ,J2,---,J n ) (2.5) 

where the z/j's are a set of constant functions of the action variables. The equation 
has the solution 

Ui = Pit + pi. (2.6) 

The constant z/j are just the frequencies associated with the periodic motion, and 
/3j's are constants of integration. This formalism ()2.6|) then gives, the frequencies of 
periodic motion. The semi-classical Bohr- Sommerf eld quantization rule is obtained 
if we require that the action- variables are integral multiples of Planck's constant. 

We will now summarize the QHJ formalism and give an exact quantization rule 
and its relationship with Schrodinger formalism. We will first give the QHJ equation 
for one dimensional system which can be easily generalized for a separable system in 
several dimensions. 

2.2 Quantum Hamilton- Jacobi Equation 

In the quantum theory, ( with 2m = 1), one assumes the generating function W(x, E) 

Sett lsflGS 

a^,fl ) + [CTr(»,fl)j' = 



i dx 



2 



dW{x,E) 



dx 

which will be called the quantum Hamilton- Jacobi (QHJ) equation. The momentum 
function 

will be called the quantum momentum function (QMF). In the limit h — > the QHJ 
equation goes over to the classical Hamilton- Jacobi equation (|2.2j) . Also the QMF 
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tends to classical momentum function 

p(x,E)^p cl (x,E) = [E-V(x)\* (2.9) 
From (|2.7jl it is seen that the QMF satisfies the following equation 

p 2 (x, E) - ihp'(x, E)-[E- V{x)\ = (2.10) 
This equation will also be referred to as the QHJ equation. 

2.3 Exact Quantization 
2.3.1 Boundary condition 

The QHJ equation ()2.1(Jj) is a non-linear differential equation and will give rise to two 
solutions. We need to establish a boundary condition to select physically acceptable 
solution. We first summarize the original boundary condition proposed by Leacock 
and Padgett. We state the boundary condition, which will complete the definition 
of the QMF p(x,E) in terms of the classical momentum function p c i(x,E). The 
classical momentum function p c i(x, E) defined by ()2.9|) is a multi- valued function of 
x and is defined by the following rule: 

The turning points X\ and x^ are defined by the vanishing of p c i(x,E) i.e by 
Pa(xi,E) = p d (x2,E) = 0. The complex x plane on which p d (x,E) is defined, is 
given a cut connecting the two branch points, i.e., a cut from x\ to x<i- p c i{x, E) is 
defined as that branch of the square root, which is positive along the bottom of the 
cut. 

With the above definition of the classical momentum function p c i{x, E), we state 
the physical condition which completes the definition of the QMF p(x, E) as: 

p(x,E) *^p d ( x ,E) (2.11) 
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Requirement (|2.11l ) has two interpretations: (1) as a form of the correspondence 
principle and (2) as a boundary condition on p(x, E). 

The boundary condition given above is easy to implement only for very simple 



potentials because p c i(x, E) = J E — V(x) will in general have several branch points, 
and the correct branch need to be selected. For this reason, we will impose other 
condition to select the solution. Several possibilities exists for an alternate condi- 
tion. For example, the square integrability of the bound state wave-function is one 
such requirement. Some other conditions, useful in the context of super symmetric 
potentials models, can also be written down [4]. 

2.3.2 Exact Quantization Condition 

Having introduced the QMF p(x, E), we define the quantum action variable by gen- 
eralizing the classical definition. The classical action variable can be defined as the 
integral 



2vr Jc 

where the integral is around a closed contour C. The contour C encloses the cut of 
Pd(x, E) which runs between the turning points x\ and x^- 

Following the above definition, we define the quantum action variable by 



where p(x, E) is the quantum momentum function, and C is the contour defined 
immediately above. 

The definition ()2.12|) connects the action-variable eigen-value J to the energy 
eigen-value E. In order to use ()2.12|) it is necessary to obtain the eigen-values J. 
Equation (j2.11j) and (|2.12j) imply that p(x, E) has poles of residue — ih on Re x-axis 





J 




(2.12) 
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between the turning points x\ and %2- For the ground state, first excited state, second 
excited state ■ • ■, p(x, E) has zero, one, two ■ • •, poles respectively in the potential 
well. The number of poles of p(x, E) in the potential well gives the excitation level 
of the system. Since these poles of p(x, E) are enclosed by the contour C, we have 

J = nh = J(E) (2.13) 

where n = 0, 1,2, • • •, and E is the energy eigen-value that is correlated with the 
values of nh for J. 

Equation ()2.13)1 can be inverted. Thus one has 

J = J(E) or E = E(J) (2.14) 

2.4 Connection with Schrodinger Equation 

To bring out an equivalence between the QHJ equation and the Schrodinger equation, 
Leacock defines the wave-function as 

ip{x, E) = exp [jrW{x, E)^j (2.15) 

The wave-function i/j(x, E) satisfies the correct Schrodinger equation and the appro- 
priate physical boundary conditions. That the wave-function satisfies the correct 
Schrodinger equation can be seen from the above definition and the QHJ equation 
for W(x,E). The wave-function for the bound states in one dimension has nodes 
whose number increases with energy; the wave-function for the n th excited state have 
n nodes in the classical region. Since the QMF is 

p(x,E) = -ih^ (2.16) 

these nodes are reflected as poles in the QMF and the residue of QMF at each pole 
is —ih. Therefore, if we take a contour integral 

fp(x,E)dx, 



CHAPTER 2. QUANTUM HAMILTON- J A COBI FORMALISM 



14 



along a contour enclosing the poles of QMF corresponding to the nodes of the wave 
function, we will have 

<j> p(x, E)dx = 2ni [sum of the residues] = nh (2-17) 

This quantization condition along with the singularities and knowledge of residues 
of QMF is sufficient for obtaining energy eigen-values for exactly solvable models. 
For other models, approximation schemes can be developed. 

2.5 Singularities of QMF 

The QHJ equation (j2.1U|) is of Riccati form. If V(x) has a singular point, in the 
complex plane, p(x, E) will also have singular point at that location. Such singular 
points are known as fixed singular points, and will be present in every solutions. On 
the other hand, other types of singular points with locations depending on the initial 
conditions, may also be present. These singular points are known as moving singular 
points. A well known theorem states that, the moving singular points of solutions 
of Riccati equation can only be poles. This pole will corresponds to a zero of the 
wave-function. Such a pole can only be a simple pole with residue — %h. In fact if we 
substitute, assuming b ^ 0, 

p{x,E)~ , b (2-18) 

[X - XqY 

in the QHJ equation 

p 2 (x, E) - ihp'(x, E)-[E- V(x)] = (2.19) 

and if the potential is not singular at X — Xq then r must be equal to one and 
b = —ill. Thus the residues at each moving pole must be — i%. This fact will be 
utilized throughout the thesis. 
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In the next section we show how to calculate the eigen-values by taking 
oscillator as an example. 

2.6 Energy Spectrum of Morse Oscillator 

The potential energy of the Morse oscillator is 



o\ 



V(x) = A + B 2 e~ 2ax - 2B(A + -)e~ ax 

2 



with the super potential 



and s — - 



W(x) — A — Be 



The quantum Hamilton- Jacobi equation is given by (h = 1 = 2m) 



p 2 (x, E) — ip'(x, E) 



a. , 



E — A — B z e~ zax + 2B(A + -)e~ ax 

2 



We effect a transformation to a new variable 



2B 



y = — e 

a 



The quantum Hamilton- Jacobi equation in the new variable is 



P (y> E ) +iayp'(y,E) 



E-A 2 -^y 2 + (A+^)ay 



where p(y) = p(x(y)). We define (p(y, E) by 



P(y,E) = iay4>(y,E) 



Therefore (J2.24|) transforms to 

2y Ay z a z y z 4 2 



Let 
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X (y,E) = <p(y,E) + ±- (2.27) 
2?/ 



Therefore (|2.26|) transforms to 

X 2 + X + + ~^- 2 [E -A*-^f + (A + ?L)yo} = (2.2.S) 
Ay z a z y z 4 2 



X has poles at y = and there are moving poles between the classical turning points 
. We assume that there are no more poles in the complex plane other than a pole of 
finite order at infinity. 

Residue at the fixed pole y = 0: For y = we define 

X= - + a + a lV + ... (2.29) 

y 

Using ()2.29|) in ()2.28|) and equating the coefficient of 4j, yields 

1 



&! = - [l±i-J|£-A2|] (2.30) 
2 a v 

The residue has two values, and the correct value is selected by imposing the condi- 
tion given below using the super potential viz., 

hm p(x, E) = i\Jlm W(x) 
In the y variable the above becomes (set 2m = 1) 

\xmp(y,E) = iW(y) 
which yields the value of b\ as in the lime^o as 

L A 1 

bi = ~ + - 

a 2 

Hence the correct sign of b\ is to choose the negative sign in ()2.3()|) and the hence 
the value of b\ is 
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fe = 5 [i-i-^-*i 



(2.31) 



Residue at y = oo: Now we determine the residue for the pole at infinity, for 
which we effect a transformation given by 

1 



y 



t 



(2.32) 



With x(t) = x(V*) P^l 1 transforms to 



X 2 (t)-t 2 f(t) + ^ 2 + 4 
4 cr 



a 



a , 



(E - A 2 )t 2 - — + (A + -)at 



(2.33) 



We assume an expansion for x(t) as 

= d + d x t + d 2 t 2 + ■■■ 

The residue of x(t) at t = is obtained from the integral 

— <f p(x, E)dx 
2n J 



(2.34) 



(2.35) 



which in the t variable yields the residue to be d\. To determine the residue d\ we 
use (12.34)1 in ()2.33|) . Therefore ()2.33|) transforms to 



1 X 

[d +d 1 t+d 2 t 2 +- ■ ■} 2 -t 2 [d 1 +2d 2 t+- ■ -} + -t 2 +^ 

4 a z 

Equating the constant term on both sides gives 

d n = ±1 



a 



a , 



(E - A 2 )t 2 - — + (A + -)at 





(2.36) 



(2.37) 



Equating the power of t on both sides gives 

1 



d\ 



2ad n 



A + 



(2.38) 
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The correct sign for do is chosen by the condition of square integrability on the wave 
function viz., 

ip(x) = exp (i J p(x,E)dx S j 

The above integral is bounded at infinity only if do = — | 

Summarizing we have the result that for large t, \ is given by (|2.H4jl where 
do = -~, dt by (|2~88jl 

The Quantization Rule and Eigen- values: 

We shall now obtain the eigen-values by enforcing the quantization rule 

J(E) = — <f p(x, E)dx = nh 

Z7T JC 

where C is a contour enclosing the part of real axis between the turning points in 
the complex x-plane. Changing the variable to y — — exp (—ax), the corresponding 
quantization condition in the y-plane becomes 

if I 

J(E) = — f (x-7r) d y = nh 

where C is the image in the y-plane of the contour C in the x-plane, but with 
anti-clockwise orientation which compensates for the negative sign coming from the 
derivative. In addition to the moving poles, p(y, E) has a fixed pole at y = 0. Let 71 
be a small circle enclosing the singular point y = 0, and T# is a circle of large radius 
R such that it encloses all the singularities of p(x, E). See fig(2.1). 
Hence 

h R = J(E) + 7 71 (2.39) 

where J 7l is the contour integral for the contour 71 enclosing the pole y = and Ir R 
is the contour integral for the contour T^.. 

We have evaluated the contour integral J 7l and its value is —b\. The value of 
J(E) is -n. 
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Rey 



For evaluating the contour integral I-p R we make a transformation of variable by y = j 
and hence the contour deforms to a new contour T r which encloses the singular point 
at y = oo or t = 0. The value of this contour integral Ip r has been evaluated and is 
—d\. 

Hence f!2.39|) transforms to 



h r = J{E) + L 



71 



(2.40) 



Hence to obtain the energy spectra of the Morse oscillator, we equate the residue 
of the fixed poles and the moving poles to those at infinity. Hence we have 



b\ — n = —di 



(2.41) 



Substituting the values of bi and d\ we get 



-(1 - i—\ \E -A 2 \)+n = 

2 V ct VI w 2ad 



a 

A + - 
2. 



(2.42) 



which on simplification gives the desired result for energy spectrum as 



E = A 2 — (A — naf 



(2.43) 



Chapter 3 



CACULATION OF 
WAVE-FUNCTION FOR ES 
MODELS 

In this chapter we apply the QHJ formalism outlined in the previous chapter, to find 
bound state wave-functions for several exactly solvable potential problems in one 
dimension. We will show that, by making use of elementary theorems in complex 
variables, the form of QMF can be determined completely and hence the bound state 
wave-functions are easily obtained. To determine the form of QMF we begin with 
the QHJ equation (k — 1 = 2m) 

p 2 {x, E) - ip'(x, E)-[E- V{x)\ = 

where p(x,E) is the QMF continued in the complex x-plane, and is related to the 
wave-function by 

The zeros of the wave- function will appear as poles in the QMF. According to the well 
known theorems about the nodes of wave-function, the n th excited state corresponds 
to n zeros on the real line, and there will be corresponding n (moving) poles in 
the QMF and the residue at each pole will be —i as has been discussed in chapter 

20 
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2. In addition to these moving poles, there are fixed poles corresponding to the 
singularities of the potential. We will make an assumption that QMF has no 
other singularities in the finite complex plane. The QMF turns out to be 
meromorphic, and to fix its form one needs to know the behavior of QMF for large x 
in the complex x-plane. This information can be easily read from the QHJ equation 
and hence the form of QMF can then be fixed completely. In the next section we show 
how this strategy works for the harmonic oscillator. In the remaining sections of this 
chapter we give the details of the calculation of the bound state wave-functions for 
harmonic oscillator, Morse oscillator, Poschl Teller, Eckart potentials and hydrogen 
atom. For these potentials a change of variable becomes necessary and we always 
try to bring the QHJ equation in the new variable to a form as the above equation. 
We also mention that several other potentials have been studied [5] and the bound 
state wave-functions in each case agree with the known results [8]. 

3.1 Harmonic Oscillator 

The potential energy of the harmonic oscillator is 

V(x) = -mu 2 x 2 (3.1) 
The quantum Hamilton- Jacobi equation is given by (h — 1 = 2m) 

p 2 {x, E) - ip'(x, E)-[E- ^uj 2 x 2 } = (3.2) 

The QMF p(x, E) has n poles corresponding to the zeros of the wave-function, 
and residue at each of these poles is —i. It can be proved that p(x, E) has no other 
poles except at infinity [1]. 

For large x 

1 

p(x,E) ~ ±-iux (3.3) 
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and we write 



p{x,E) » ±-iux + Q(x) (3.4) 



where Q(x) is to be determined. 

The sign of p(x, E) is determined by the condition of square integrability of the 
wave-function. 

The wave-function is expressed as 

ip(x) = exp (^i J p(x,E)dx^j (3.5) 

When the above value of p(x, E) is substituted in the equation of wave-function, the 
wave-function is bounded at at large x if we choose the positive sign of hiuix. 
Hence we write the quantum momentum function p(x, E) as 

n —i 1 

p(x, E) = ^2 1 — iwx + <p(x) (3.6) 

k=l ^ "^b ^ 

where x\,X2, ■ ■ ■ ,x n are the location of n poles on the x-axis and <f>(x) is analytic 
every where and bounded at infinity. Therefore Liouville's theorem tells us, it has 
to be a constant. Hence let <b( constant. Hence the above equation becomes 

- — i 1 

p(x,E) = ^2 h -iu)x + c (3.7) 

k=l 



Substituting ()3.7|) in ()3.2|) we have 

A -i 1. 1 
y j 1 — tux + c 

k=l ^ 



Equating the power of x, we have c = 0. Equating the constant term to zero on both 
sides in equation ()3.8|) gives 

E = {n+ l -)uj (3.9) 

which is the well known expression for energy of the harmonic oscillator in our 
notation (2m = 1) 
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The sum of moving pole terms 



E 



k=l X Xk 

can be expressed as where P(x) is the polynomial 



p( x ) - n ( ;/ x k)- 



k=l 



The QMF ()3.7|) can be expressed as 



/ ^ - p '( x ) !• 



Using (J3.ll)) in (J3.2j) and on simplification yields 



(3.10) 



(3.11) 



P"{x) - uxP'(x) + nuP(x) = 



(3.12) 



We effect a transformation £ = ax where a 2 = ^. Hence equation 1)3.12)1 changes to 



P"(O-2£P / (O + 2nP(O = 



(3.13) 



The above equation resembles the well known Hermite differential equation. Hence 
on comparison, we get 

P(0 = H n {ax) (3.14) 

and is the Hermite polynomial. 

The wave-function is expressed as 



i(}(x) = exp 
and hence we have 



i J p(x, E)dx 



exp 



P'(x) .1 
' / 1— i — , + i—i=ujx)dx 



P(x) y/2 



ip(x) = H n (ax)exp ^—-ux 2 ^J . 



(3.15) 



This is the desired wave-function for the harmonic oscillator. 
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3.2 Morse Oscillator 



The potential energy of the Morse oscillator is 



o\ 



V{x) = A 2 + B 2 e~ 2ax - 2B{A + -)< ' 



with the super potential 



and 



W(x) — A — Be' 



A 

s = — . 

a 



The quantum Hamilton- Jacobi equation is given by (K — 1 = 2m) 



p 2 (x, E) — ip'(x, E) 



E- A 2 - B 2 e~ 2ax + 2B(A + ^)e~ ax 



We effect a transformation to the variable 



2B 



y = — e 

a 



The quantum Hamilton- Jacobi equation in the new variable is 

P 2 (y,E) +iayp'(y,E) 



a 2 n , , a. 



We define 



Then (1.3.211) transforms to 



E - A 2 - -y 2 + (A + ^)ay 



P(V,E) = iay<f)(y,E) 



1 ^2 , jli 1 1 it? ,2 y 2 ® 2 



2y 



4 



+ -) 2 + , -7^ + ^^-^ 2 -^ + (^ + f)H = o 

/±yZ a ZyZ I 2 



Let 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



X (y,E) = ( j ) (y,E) + — 



(3.24) 
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Therefore ()3.23|) transforms to 

X 2 + ^~ 2 +x' + ~^- 2 [E -A*-^f + (A + |)yn] = (3.25) 
Ay z a z y z 4 2 



This equation suggests that \ has a pole at y — 0. It will also have n moving 
poles corresponding to the nodes of the wave-function. We assume that there are no 
other poles in the finite complex plane. For large y the behavior of \ has already 
been obtained in section (2.6) and is seen to be bounded for large y. Hence we get 
using Liouville's theorem 

X = - + ±(—)+c (3.26) 

where b\ and c are constants to be fixed. The residue of x a t V = is b\ and has 
been obtained in section (2.6) and 

&! = - + ! (3.27) 
a 2 



We write, once again, 



where 



y i P\y) 



P{V) = U(V - Vh) 



k=l 



Substituting lEPEI) in (l3~^oT) gives 



— + 2—— + 2— c + 2— c + c 2 - - + — (A + — ) = (3.28) 
P y P P y 4 ay K T v ' 

In order to proceed further we look at the behavior of each term for large y. Using 
the leading terms 

P"(y) n(n - 1) P'(y) n 



P(y) y* ' P(y) y 

in equation (|3.28j) and equating the constant term on both sides gives, c = ±|. 

The correct sign for c is chosen by the condition of square integrability on the wave 



function which fixes c = — \. 
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Comparing the coefficient of ^ for large y on both sides we get 

26ic + 2nc + (a + -) - = 
\ 2/ a 



(3.29) 



which on using the values of b\ and c and on simplification gives the energy eigen- 
value 

E = A 2 - (A - naf (3.30) 



Substituting the value of &iand c in equation (|3.28|) we have 

yP'\y) + {l-y + 2(s- n)}P'{y) + nP(y) = 
Compare this with the standard Laguerre differential equation 

xy" + (P + 1 — x)y' +ny — 

we have P(y) = L^(y). 

The wave-function for the Morse oscillator is given by 



ip(x) = exp J p(x,E)dx 



In terms of y variable we have 

ip(y) = exp (^i 
On integrating and simplifying we get 



bt P'(x) 1 1 
7 + T{x) ~ 2 " Yy 



dy 



(3.31) 



(3.32) 



(3.33) 



My) =y s "ex P (--?/)P(x) 



Replacing the value of P(y) we have 



(3.34) 



My)=y s ~ n zM-~y)L n p{y) 



(3.35) 
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3.3 Poschl- Teller Potential 



The Poschl- Teller potential is 



V(x) = A 1 + (B 2 +A+ Aa) csdraz - B (2A + a) coth ax cschax (3.36) 



with the super potential given by 



W(x) = A coth ax — B eschar (A < B) 



(3.37) 



and 



A (3 
s = — , A = — 



(3.38) 



a a 

The quantum Hamilton- Jacobi equation is given by (h = 1 = 2m) 

p 2 (x, E) — ip'(x, E) 

- E - A 2 - (B 2 + A 2 + Aa) csch 2 ax + B(2A + a) coth ax csch ax = (3.39) 
We effect a transformation to a new variable 



y = cosh ax 



(3.40) 



The quantum Hamilton- Jacobi equation in the new variable becomes 



p (y, E ) - ia yy 2 - W(y, E) 

E-A 2 -(B 2 + A 2 + Aa)^— + B(2A + a" U 



y 2 -i 



y 2 -i 



We define 



(3.41) 



P(y,E) = -ia^Jy 2 - l<f>(y). 
Therefore (13.4111 transforms to 



(3.42) 
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+ 



Let 



a 2 (y 2 — I) 



E-A 2 -{B 2 + A 2 + Aa)^— + B{2A + a" U 



y 2 -l 



y 2 -l 



X = 4> + t;-^— 



2y 2 



(3.43) 



(3.44) 



Therefore ()3.43|) transforms to 

x 2 + x' + 



3 y 2 



1 1 



4(y 2 



+ 



1 



a 2 (y 2 — 1) 



E- A 2 - {B 2 + A 2 + Aa)- 



1 



2y 2 



+ B(2A + a) 



y 



(3.45) 



v - 1 v 2 -i_ 

X has poles at y — ±1 and there are moving poles between the classical turning 
points. We assume that there are no more poles in the complex line. We now 
determine the residue at each of these poles. 
For y = +1 , we define 

6i 



X 



+ a + ai(y - 1) + 



(3.46) 



Using (|3.46|) in (|3.45|) and equating the coefficient of 7~jp yields 



1±— [2(S-A)-a] 



(3.47) 



The correct value of &i is selected by imposing the condition 



lim Res p(y, E) = iRes W(y, E) 

E — >0 



on super potential as explicitly shown for Morse oscillator gives 



(3.48) 



For y 



T , we define 



X 



y+1 



+ a' Q + a[(y + I) + 



(3.49) 
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Using ()3.49|) in f)3.45|) and equating the coefficient of t^jw and following the above 
procedure yields 



(3.50) 



The residue at a moving pole is seen from ()3.45j) to be 1. Hence we arrive at the 
form 



X 



h b\ , P'{y) 



+ c 



y-1 y + 1 P(y) 
for x where c is a constant to be determined 

Substituting (|3.51|) in (|3.45j) gives the following equation. 

, 26i&'i 2b\ P' P' 2b lC 2b x P' 2b\ 1 1 P" 

c + ^T^T + V^7 + 2 ^7 C + 7 + 7^7 + ~T C ~ n~? 7 + "77 

y 2 - 1 y + lP P y-1 y - 1 P y + 1 2y 2 -1 P 

3 1 



(3.51) 



1 :(B 2 + A 2 + Aa)- ' 







(3.52) 



8y 2 — 1 a 2 y 2 — 1 2a 2 y 2 — 1 

Now we look at different terms in equation ()3.52j) for large y and equate their coeffi- 
cient to zero. Equating the constant term to zero gives c = . With c = the above 
equation becomes 

2b 1 b[ 1 1 



3 1 



■ y 2 — 1 a 



P" P' 

~p + T 



26i 2b\ 

— H — 

y-1 y+l 

, 1 1 



y 2 - 1 2cr 



y 2 - 1 2y 2 -l 



(B 2 + A 2 + Aa)- 







(3.53) 



P"(V) _ n(n-l) _ n 



For large y, P(y) behaves as P(y) ~ y n H and -p^ ^— , p(j/) - y 

Using these in f!3.53|) and equating the coefficient of p gives 

2M'i + 2nb[ + 2nh + n(n - 1) + - + — (£ - A 2 ) + -^(£ 2 + A 2 + Aa) - - = 



2c^ 



and substituting the values of b\ and b' x gives the expression for energy as 
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E = A 2 — (A — no) 2 (3.54) 
Substituting the value of E, b\ and b[ in ()3.53|) we get 

(1 - y 2 )P"(y) + [(2s - l)y - 2\]P'(y) + [n(n - 2s)]P(y) = (3.55) 

The above equation resembles the standard Jacobi polynomial, such that 

P(y) = P<*P>(y) = Pi X - s -^ X - s ^\y) (3.56) 

The wave function for the Poschl Teller potential is obtained on the same line as 
that for the Morse oscillator and is given by 



^( y ) = ( y -l) L ^ y + l)^pi x 2 ' 2 \y) 



which agrees well with the values given in literature. [5] 



(3.57) 



3.4 Eckart Potential 



The potential energy of the Eckart potential is 

B 2 



V(x) = A + — — 2B coth ax + A(A - a) csclTax 
A 1 



with the super potential given by 



B 



W{x) = -Acothax+ — , (B > A 2 ) 



and 



A 

s = — 

a 



a=4 

a 2 



A 



a 



n + s 



The quantum Hamilton- Jacobi equation is given by (H — 1 = 2m) 

B 2 



p (x, E)—ip'(x, E) — 



E — A — — + 2B coth ax - A(A - a) csclTax 



(3.58) 



(3.59) 



(3.60) 



(3.61) 
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We effect a transformation by the variable 



y = coth ax 



(3.62) 



The quantum Hamilton- Jacobi equation in the new variable is 

B 2 



p 2 (y,E)-za(l-y 2 )p'(y,E) 



We define 



E - A 2 - — + 2By - A(A - a)(y 2 - 1) 



P(y,E) = -ia(l - y 2 )<p{y) 



Hence equation ()3.63|) simplifies to 



y i2 



a 2 (l — y 



2\2 



l-2/ 2J 
B 2 



+ 0' 



J--y 



2\2 



Let 



E - A 2 - — + 2By - A(A - a)(y 2 - I) 



x = <P- 



Therefore the above equation changes to 



(3.63) 



(3.64) 



(3.65) 



(3.66) 



+ 



1 



X +X + 



B 2 



1 — y 2 ) 2 I — y 



E - A 2 - — + 2By - A(A - a) (y 2 - 1) 







(3.67) 



a 2 (l — y 2 ) 2 

X has poles at y — ±1 and there are moving poles between the classical turning 
points. We assume that there are no more poles in the complex plane. We determine 
the residue at each of these poles. 
For y = +1, we define 

h 



X 



+ a + a x {y - 1) + • 



(3.68) 
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Using ()3.68|) in ()3.67|) and equating the coefficient of ( y _\)2 yields 



61 = - 
1 2 



(3.69) 



As the residue has two values, the correct value is selected by imposing the condition 
on the super potential as done in Morse oscillator and the correct value is 

1 



1 + 



a 



(3.70) 



Similarly the residue for y — — 1 is determined and is given as 



We assume x t° have the form 



X 



' + — + -TTT^f + c 



(3.71) 



(3.72) 



2/ - 1 y + 1 P(y) 

where c is a constant to be determined 

Using ()3.72|) in ()3.67|) and following the similar lines as that of Morse and Poschl 
Teller potential one gets the value of c = and the resulting equation becomes 



1 1 



P"{y) P'(y) 
P(y) P{y) 

1 1 



261 b[ 



y - 1 y + l 



+ 



2M'i 

y^i 



\E-A 2 -%-\- ^-^—A(A -a) = (3.73) 



2y 2 -l 2a 2 y 2 -V A 2> a 2 y 2 - 1 

For large y assuming P(y) ~ y n + • ■ ■ and equating the coefficient of 4? gives 

26 1 6' 1 + 26> + 2M + n{n - 1) - J - ^{E -A 2 -%-)~ \a{A - a) = 

2 2ar /r cr 

Substituting the values of b\ and 6^ gives the energy expression as 



E = A 2 - (A + naf - 



B 2 



B 2 



(A + na) 2 A 2 



(3.74) 
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Using the values of b±, b[ and E in ()3.73|) one gets the differential equation for Eckart 
potential as 

(1 - y 2 )P"(y) + [-2a - 2(-n - s + l)y] P'(y) + (-2ns)P(y) = (3.75) 
Equation ()3.75j) resembles the standard Jacobi polynomial and 

P{y) = P^\y) = P^\y) (3.76) 
The wave function for the Eckart Potential is obtained from 

and is given by 

1>(y) = {y - 1)% (y + l)^P^\y) (3.77) 
The values for energy and wave function agree with those found in the literature. [5] 

3.5 Hydrogen Atom 

In this section we obtain the bound state wave functions of the radial part of the 
Schrodinger equation (H = 2m = 1) 

d 2 R | 2 dR | / | Ze 2 A \ R _ Q 
dr 2 r dr \ r r 2 J 

where A = /(/+1). Using the transformation R(r) = 0(r)/r, the Schrodinger equation 
becomes 

^ + (£+ ££!_A W = . (3.79) 



dr 2 r r 



The QHJ equation in terms of 



g=^ln(0(r)) (3.80) 



CHAPTER 3. CACULATION OF WAVE-FUNCTION FOR ES MODELS 



34 



is given by 

92 + ^ + ( B + ^!_A) =0 . (3.81) 
ar \ r r A ) 

The range of r is from to oo, and the wave function <f>(r) should vanish at r = 0. 
Thus q has a fixed pole at r = 0, along with the n moving poles with residue equal 
to one on the real line. Like harmonic oscillator there are no other singularities in 
the finite complex plane. Thus we can write q, in a similar fashion as for harmonic 
oscillator, as 

q(r) = L + h + C (3.82) 
P r 

where b\ is the residue at r = which can be obtained by doing a Laurent expansion 
of q around the pole at the origin. The two values of residues obtained are 

b x = -l, b x = l + 1. (3.83) 

One chooses the right residue by using the square integrability property of the wave 
function (p and obtain 

6i = I + 1 (3.84) 

as the right choice. Substituting (|3.82j) for q in ()3.8H) and expanding different terms 
of the resulting equation for large r and comparing the leading terms we get 

C 2 = —E, E = --—-^ (3.85) 
(2n') 2 K ' 

where n' — n + I + 1 and one is left with the differential equation 

r r + 2P ( l + 1 .^L) + ("'-'- 1>^p = , (3.86) 



Now defining 



2n' I n' 



7 ^ r r = p (3.87) 

n 



()3.85|) becomes 

pP" + ((2/ + !) + !- p)P' + (n' - I - l)P = (3.88) 
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which is the associated Laguerre differential equation where P is the Laguerre poly- 
nomial denoted by L. The bound state wave function obtained from (J3.82j) and 
(jggUj) is 

MP) = P l+1 exp(-p/2)L»+}(p) (3.89) 
which is seen to be identical with known correct answer. 



Chapter 4 



CONDITIONS FOR 
QUASI-EXACT SOLVABILITY 

4.1 Introduction to QES 

In this chapter we study QES model in one dimension. These are the models for which 
a part of the bound state energy spectrum and corresponding wave-functions can be 
obtained exactly. These models have been constructed and studied extensively by 
means of Lie algebraic approach. For a review we refer to the book by Ushveridze et al 
[6]. In order that a part of the spectrum be obtained exactly, the potential parameters 
appearing in the potential must satisfy a condition known as the condition for quasi- 
exact solvability. Within the QHJ approach, as used for exactly solvable models, it 
is not clear how such a condition can arise and why only a part of the spectrum is 
exactly solvable. In this chapter we report a study of these aspects of QES models. 

In order to study QES models within QHJ formalism one needs to have informa- 
tion of singularities of QMF. This in general is not very easy to obtain except for 
some simple cases like harmonic oscillator and hydrogen atom problems. In the limit 
h — > the QMF p(x, E) goes over to p c i(x, E) = \J E — V(x) which will, in general, 
have several branch points. This is an indication that in general, the singularity 
structure of p(x, E) will be very complicated. In order to make progress, we 

36 
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make a simplifying assumption that the point at infinity is an isolated 
singular point and more specifically it is a pole of some finite order. Thus 
this amounts to saying that p(x, E) has fixed poles, and a finite number of moving 
poles and a pole at infinity. Using these requirements, we can proceed as in the 
case of exactly solvable models and work out the consequences of exact quantization 
condition given below. 



We find that for all the QES potential models studied by us, (j4.1|) and our 
assumptions, imply that potential parameter must satisfy a condition which turns 
out to be identical with the condition of quasi-exact solvability of the potential. A 
list of potentials studied and the condition of quasi-exact solvability in each case are 
given below. 

The potentials are: 

1. Sextic oscillator: 




(4.1) 




7 > 0. 



(4.2) 



2. Sextic oscillator with centrifugal barrier: 



V{x)=4(8-±)(s 




3. Circular potential: 




sin x 



COS z X 



(4.4) 



with 



A = 4{ 8l -\)( 8l -l), 



(4.5) 




(4.6) 
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C = q\ + Aq 1 (s 1 + s 2 + fi), (4.7) 

D = q\. (4.8) 

4. Hyperbolic potential: 

A B 

V(x) = + — n C cosh 2 x + D cosh 4 x, (4.9) 

cosh x sinh x 

with 

A = 4( ai -l)( ai -^), (4.10) 

B = 4( 82 -\)(82-l), (4.11) 

C = [g 2 + 4 gi ( Sl + S2 + /i)], (4.12) 

D = q\. (4.13) 

5. 

V(x) = A sinh 2 y/ux + B sinh y^x + Ctanh v / z/xsechy / z/a; + Dsech 2 v / E'a; (4.14) 

6. 

V(x) = A cosh 2 + B cosh + C coth yfvx esc hy/vx + Dcsch 2 A/z/a; 

(4.15) 

7. 

= Ae 2 ^ x + Be^ x + Ce"^ 7 * + De' 2 ^ x (4.16) 
The conditions for quasi exact solvability for these potentials are: 
1. 

1 ( I3 2 \ 

a] = 3 + 2n, n = integer (4-17) 

y/l V 4 7 / 
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2. 



/j, = integer 



(4.18) 



3. Taking 



A = A(s 1 - 1 -)(s 1 - 3 -) 
B = A{s 2 -\){s 2 -\) 
C = ql + 4 qi ( Sl + s 2 + /i) 
D = q\ 



(4.19) 
(4.20) 
(4.21) 
(4.22) 



the condition is 



pi = integer 



(4.23) 



4. Taking 



A = 4( ai -i)( ai -|), 
B = A{s 2 -\){s 2 -\l 
C = [qf + A qi ( Sl + s 2 + //)], 

£> = ?i 2 , 



(4.24) 
(4.25) 
(4.26) 
(4.27) 



the condition is 



/i = integer 



(4.28) 



5. 



B±2(n + l)VuA +A(4D-u) B ± 2(n + 1) yJvA - 4A 2 C= (4.29) 



6. 



B±2(n + l)VvA -A(4D + v) B±2(n+l)VvA + AA 2 C 2 = (4.30) 
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7. 

2(n + l)y/vAD = ±bV15 ± Cv 7 ^ (4.31) 
The calculations for these potentials are given in the next few sections. 

4.2 A Representation of QES Quantization Rule 

We now bring out some common features of the exactly solvable models, that have 
been studied in this thesis and those reported in the paper [5]. For the exactly 
solvable model the QMF written in terms of suitable variables y takes the form 

p{y) = + -\ + ••• + £ — L- + R(y) (4.32) 

where £1,^2, • • • are fixed poles, the summation term corresponds to n moving poles 
at 2/1, 2/2, • • • , Vn and R(y) is atmost a polynomial in y. The residues 61, b[, • • ■ have 
been calculated using the QHJ equation and demanding a condition such as the one 
proposed by Leacock and Padgett, or 

lim p(x, E) = i W(x), 

E — >0 

or the square integrability of the wave-function. Thus in all the cases studied we are 
lead to a rational expression for the QMF. 

Under an assumption about the behavior of QMF at infinity, even for QES mod- 
els, the QMF turns out to be a rational function. The quantization rule 

£p(x,E)dx = nh (4.33) 

is then easily seen to be equivalent to the well know result, that for a rational 
function, sum of residues at all poles, including the one at infinity vanishes. Written 
explicitly for a rational form of QMF that we have, this requirement becomes 

(Res of QMF) + n + (Res of QMF at infinity) = (4.34) 

fixed poles 
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where Res stand for the residue and the middle term n, corresponds to the contri- 
bution of moving poles to the residue. 

In this and the next chapter, we will use this condition (|4.34|) as a substitute for 
quantization rule. 

4.3 Sextic Oscillator 

The potential for the sextic oscillator is: 

V(x) = ax 2 + px 4 + jx 6 , 7>0 (4.35) 

The QHJ equation is (h = 1 = 2m) 

p 2 (x, E) - ip'(x, E) - (E - V) = (4.36) 

For the n th excited state, the QMF has n poles on the real axis and we assume 
that there are no other moving poles. We shall use the quantization condition viz., 

— / p(x, E)dx = nh (4.37) 

Z7T JC 

in the form (|4.34|) as given above. 

To evaluate the integral in ()4.37|) a Laurent expansion of p(y) in powers of y = 
1/x, is made 

P(V) = + - + a o + a lV + ■■ ■ (4.38) 

yi yZ y 

Substituting this in (|4.36J) and integrating term by term we get 

J(E) = tax (4.39) 

The quantization condition gives 

a x = -in (4.40) 
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It only remains to compute the coefficient a\ of the Laurent expansion given in (|4.38j) . 
To do this we start from the QHJ equation 

p 2 (y) + iy 2 p'{y) -E + ^ + ^ + ^ = (4.41) 

Substituting the Laurent expansion and equating the coefficients of different powers 
of y on both sides of the equation we get 

h = ±i^/l ( 4 - 42 ) 

b x = -^b 3 (4.43) 

ai = + (4 . 44) 

2o 3 

It is important to know that, we would get two solutions for b\ corresponding 
to the two solutions of b 3 = ±2^/7. This happens due to the fact that the QHJ is 
quadratic in the QMF. Thus one needs a boundary condition to pick the correct 
solution. We propose to use the square integrability of the wave-function instead 
of the original boundary condition, explained in chapter 2, which was proposed by 
Leacock and Padgett. This is because the original boundary condition is difficult to 
implement in the present case due to the presence of six branch points in the p c i. In 
order to find the restrictions coming from the square integrability, we compute the 
wave-function 

ift(x) = exp ^ ip(x)dx^j (4.45) 

for large x as follows. The most important term in the Laurent expansion (14.38)1 for 
small y ~ 0, corresponding to large x is 

P(y) ~ (4-46) 
and the wave-function for large x becomes 

ip(x) w exp (4.47) 
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Out of the two solutions, 63 = ±2^/7, if)(x) is square integrable only for 63 = iy/j. 
Using this value of 63 and from (J4.43)) and (|4.44jl . equating a\ to —in we get 



a ) = 3 + 2n (4.48) 



1 ((3 2 

I 

In order to compare the results in 1)4.48)1 with the well known condition, we write 

7 = a 2 , p = 2ab (4.49) 
Thus we get a = b 2 — a(3 + 2n) which agree with the result given in [6]. 

4.4 Sextic Oscillator with a Centrifugal Barrier 

The potential is given as 



V(x) = 4(s - j)(s - ^)^2 + [b 2 - 4a(s + ^ + fx)}x 2 + 2abx 4 + a 2 x 6 (4.50) 



We shall consider only the case s > | so that the coefficient of the centrifugal term, 
^2 is positive. The Q.H.J equation is (h = 1 = 2m) 

p 2 (x,E) -ip'(x,E) - (E- V) =0 (4.51) 

Substituting the potential the QHJ equation is 

p 2 {x,E)-ip'{x,E)-\E-4{s--){s-h^--[b 2 -Aa{s + - + fi)]x 2 -2abx 4 -a 2 x 6 = 

4 4 x 2 2 

(4.52) 

p(x, E) has poles at x = and as the potential is symmetric there are moving poles 
on either side of the origin. We assume that there are no more poles in the complex 
plane. We assume that infinity is a pole. We find below the residues for each of these 
pole. 

We expand p(x, E) as 

p(x, E) = — + ao + ajx + • • • (4.53) 

x 
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Using ()4.53|) in (|4.52|) and equating the coefficient of \ we get 

h = -- [1 ± (4a - 2)] (4.54) 
Demanding that the wave-function remains finite, for x — > 0, gives 



61 = ~\ [4* - 1)] (4-55) 

To find residue for the pole at infinity, we use the mapping x — |. Therefore 
equation ()4.52|) transforms to 



p 2 {t,E)+it 2 p'(t,E) 

1 1 _ 1 i 

(4.56) 

± z i»~ v u" j 

We expand p(t, £") as 



£ - 4( a - - ^)i 2 - [6 2 - 4a( S + \ + „)]! - 2a&l - a 2 l 



(La o _n 

p(t, £) = -± + -± + -j + c + cit + c 2 t 2 + ■■■ (4.57) 
Using ()4.57|) in ()4.56|) and equating different coefficient of t we get 

d 3 = ±ia (4.58) 



—ab 

d 2 = -r- (4.59) 
c = (4.60) 

The correct sign of d 3 is fixed by the condition of square integrability and is given 
by 

d 3 = -ia (4.62) 
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Now equating the sum of all residues to zero, we get 



ib\ + 2n = ci 



(4.63) 



Substituting the values of &iand C\ in the above relation yields the required condition, 
viz 

n = fi (4.64) 
The above condition agrees with those given in [6] 

4.5 Circular Potential 



The potential is given as 

V{x) 

where 



A B . Q . 4 

H - hCsm x — JJ sin a; 



sin 2 x ' cos 2 x 



A = A(s 1 - 1 -)(s 1 -- 4 ) 
B = A(s 2 - 1 -)(s 2 - 3 -) 
C = q\ + 4gi(si + s 2 + /i) 



= 



(4.65) 

(4.66) 
(4.67) 
(4.68) 
(4.69) 



We effect a change of variable by 



y — sin x 



(4.70) 



The Q.H.J equation is 



p 2 (x, E) - ip'(x, E) - (E - V) = 



(4.71) 



In the new variable the QHJ equation is 



p 2 (y,E)-2iy/yJl-ytf(y,E)- 



A B , 
E Cy + Dy 2 

y i -y 



(4.72) 
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Let 



p(y,E) = -ziVyJ 1 - v<f> 



(4.73) 



In terms of the above transformation the QHJ becomes 

A B 



., ,, l(l-2y), 1 1 



2y(l-yY 4^/(1-1/) 



(4.74) 



y 1-2/ 

has poles at y = and at y = +1 and there are a finite number of moving poles in 
the complex plane. We assume that there are no more poles in the complex plane. 
We find the residues for each of these pole below. 
For y = we consider an expansion in as 



y 



a + aiy + 



(4.75) 



Using ()4.7f)j) in 1)4.74)1 and equating the various powers of y we get the following. The 



power of \ gives 



\ ± (2ai " 1) 



(4.76) 



The correct value of b\ is fixed by the condition of square integrability of the wave- 
function and is given below as 



+ (2si - 1) 



(4.77) 



For y = 1 we consider an expansion in <ft as 
, ( b[ 



JJ 



- + 0o + ai(2/- 1) + 



(4.78) 



Using ()4.78)1 in ()4.74|) and equating the various powers of y we get the following. The 
power of gives 



± (2s 2 - 1) 



(4.79) 
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The correct value of b[ is fixed by the condition of square integrability of the 
wave-function and is given below as 



1 2 



\ + (2* 2 - 1; 



(4.80) 



To find residue for the pole at infinity, we use the mapping y — -. Therefore 
equation (|4.74|) transforms to 

2 (t, E) - t 2 4>'(t, E) + E ) 



E t 2 At 3 



2 t- 1 
B t 3 C t D 1 



At- I At- I 4(t-l) 2 At- I At-1 
We expand (fi(t,E) as 

4>(t, E) = ( j + c + c 1 t+--- 



(4.81) 



(4.82) 



Using ()4.82|) in ()4.81|) and equating the power of i we get 



di = Q (4.83) 

Equating the constant term we have 

co = ±| (4.84) 
Equating the coefficient of t gives the residue at y = oo as 

Cl = 2^ 2 (4 ' 85) 

The correct sign of Cq is fixed by the condition of square integrability and is given 

by 

co = | (4.86) 
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Now equating the sum of all residues of fixed poles and the moving poles and the 
pole at infinity, we have the following relation. 

&i + b[ + n = ci (4.87) 

Substituting the values of bi,b[, and c\ in the above relation yields the required 
condition 

H = n (4.88) 
The above condition agrees with those given in [6] 

4.6 Hyperbolic Potential 

The hyperbolic potential is 

A B 

V(x) = + — n C cosh 2 x + D cosh 4 x (4.89) 

cosh x sinh x 

where 

A = 4( ai -l)( Sl -|) 
B = 4(s 2 - 1 -)(s 2 - 3 -) 
C= [g 1 2 + 4g 1 (s 1 + s 2 + / u)] 
D = q\ 

We will consider the case s 2 > f ■ 
The Q.H.J equation is (/i = 1 = 2m) 

p 2 (x,E) - ip'(x,E) - (£-F) =0 (4.90) 

We use a mapping is y = cosh x 

The Q.H.J equation in the new variable is: 
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p(y, E) - ih^fy* - lp'(y, E) - [E + - - + Cy 2 - Dy 4 ] = 

y y i 



Let 



p(y,E) = -iyJy2-l(l>(y,E) (4.91) 

and hence 



V \y, E) = _ 1(j) > + 0^#=} (4.92) 



Therefore the Q.H.J, equation becomes 

2 + 0' + + -J-r[E + 4 " "T^T + ^ " D V*\ = ( 1035 

y — l y — i y y — 1 



We note that has fixed poles at y = ±1 and y = and we write the Q.H.J 
equation as 

, 2 ,/ rl 1 1 1 , , . „rl 1 11, . r 1 11 11 



+ ^ + to — rr + o — + E h — t ~ n—r) + A l~^ - + 



2y+l 2y-Y T 2y-l 2y + V 1 y 2 2y + l 2y-Y 

„r 1 1 11 11 11, 

- fl -7— -r + 77—^ + + 



4y-l 4(y-l) 2 Ay + l 4(y + l) 2J 

+ C1+ -Dl + y 2 + =0 4.94 

L 2y-1 2y+l J 1 y 2y-l 2y+l J V ; 

For y — 
Let 

= - + E a ^ n 

2/ n=0 



~i / ^ = '2^f {b j + t^y n )dy = -^2nih = h 

The QHJ equation (|4.94j) becomes 

,b\ , 6i r 1 1 1 1 , , „ r l 1 11, 

(- + a + a 1 y + ...) 2 + (-- + a 1 + ... ) + [-—— + - }(/) + E[- - -— — 

V V 2 2y + l 2y-l 2y-l 2y + l 
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111 1111 11 11 11 

+ ^~f~2y + l + 2y-V [ ~4 y - 1 + 4 (y - l) 2 + l~y + l + A{y+l)^ 

11 11, r 9 11 11, , 

+ C1 + -Dl + )/ 2 + =0 4.95 

1 2y-l 2t/ + l J [ 2y-l 2t/ + l J 1 ; 

The coefficient of ^ gives 

b\ - bx - A = (4.96) 



and hence 



For y = 1 , let 



Therefore 



&i = \ ± (2si - 1) (4.97) 



y oo 

^T + E^d/- 1 )™ (4-98) 



n=0 



CO 



Substituting ()4.98j) in ()4.94j) and equating the coefficient of jzhm we get 



which yields 



,/2 1., -1 



For y 



-1, let 



6i = - ± -(2s a - 1) 



(4.100) 



ft" 00 

-XT + E<(2/ + l) n (4.101) 

y T 1 n 



n=0 



Therefore 
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Substituting (|4.101J) in (|4.94|) and equating the coefficient of r-fW we get 



(y+l) 2 

bf - -b'l - - = (4.102) 



On simplification yields the same value as b[. Hence 



&? = 1±~(2 S2 -1) (4.103) 

To calculate the residue at y = oo, we apply the mapping y = - 
The QHJ equation ()4.94|) becomes 

~ ~ I II 1 II ~ 

<f> 2 (u, E) + (-u^'iu, E) + [-—+ 273H^> E ) 
+ E[^^- 1 -^}+A[-u 2 1 U ■ 1 U i 



L 21 + u 21 + u J 21 + m 21-u 

1 U 1 V 2 1 It 1 M 2 . 1 U 1 U 

~ B [-1-, + 777—77 + 77— + 777-7-7?] + C \ l + 



41 — u 4(l-u)2 41 +m 4(1 + m) 2J l 2 1 -u 21 + w 

11m In 
+ 2 1 -u ~ 21 + W 



-£>[! + _ + -- _____] = o (4.104) 



Let 



(u,E) = — + d + d 1 u + --- (4.105) 

u 



The residue at infinity is the coefficient of d\ and is given as 

d, = - 2gl(Sl + S2 + /i) - 1 (4.106) 

Cl 

with 

Cl = ± 9l (4.107) 

Equating the sum of all residues of fixed poles, the moving poles to the pole at 
infinity, we have the following relation. 

d 1 = 2n + b 1 + b' 1 + b'i (4.108) 
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The equation (4.6.6) does not change when replacements y — > — y and — > —0 is 
made. Therefore we select the residue b[ = b'{ and the condition of finiteness of the 
wave-function at x = requires that positive sign be selected in (4.6.12) and (4.6.15). 
Therefore 



The point y — corresponds to complex value of x and therefore one cannot insist 
on finiteness of the wave-function at x — 0. One must fall back on the boundary 
condition given in chapter 2 section 2.3.1. We will simply note that selecting positive 
sign in (4.4.19) leads us to the correct condition 



for quasi exact solvability [6]. In this chapter, our objective has been to show that 
QES conditions follows from our assumption that point at inifnity is an isolated 
singular point. In the cases where for certain ranges of potential parameters, both 
the residues are acceptable, one must accept all such answers and work out the 
consequences. This may lead to some new and interesting results as is evidenced by 
the investigations on phases of super symmetry [13] and periodic potentials [15]. 

Besides the above QES potentials, we now take up three classes of QES potentials 
[9,10] and find the conditions for quasi exact solvability within our approach. 



b[ = K = S2~\ 



(4.109) 



ji — n 



(4.110) 



4.7 V(x) = A sinh 2 y/vx+B smb. ^JDx+C tanh y/vx sechy / z/x+ 




The potential is given as 



V(x) = A sinh 2 \fvx + B sinh \J~ux + C tanh ^/uxsech^/ux + D sec h?\/vx (4.111) 
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We effect a change of variable by 



y = sinh \pu. 



vx 



(4.112) 



The QHJ equation is (h — 1 = 2m) 



p 2 (x, E) - ip'(x, E) - (E - V) = 



(4.113) 



In the new variable the QHJ equation is 



p (y, E ) - iy/vy i + y 2 P'(y, E) 



.2 n.. r< y _ D _ 1 



E — Ay — By — C 



1 + y 2 1+y 2 



Let 



P{y,E) = -iy/uJl+y 2 



= 
(4.114) 

(4.115) 



In terms of the above transformation the QHJ becomes 



[0 + 



i y 



21 + y 2 



+ 0'- 



i y 2 



4 (1+y 2 ) 2 



+ 



v(l+y 2 ) 



y 



E — Ay — By — C-f — - D 



1 



Let 



X = + 



1 + y 2 1 + y 2 
2 1 +y 2 



(4.116) 



(4.117) 



Therefore the above equation becomes 



X + X + 



3 y 2 



1 1 



4 (1 + y 2 ) 2 2 1 +y 2 



+ 



v(l+y 2 ) 



E — Ay 2 — By — C- 



y 



D- 



1 







(4.118) 



1 + y 2 l + y 2 _ 

X has poles at y = ±i and there are moving poles between the turning points. We 
assume that there are no more poles in the complex line. We find the residues for 
each of these pole below. 



CHAPTER 4. CONDITIONS FOR QUASI-EXACT SOLVABILITY 



54 



For?/ = % we consider an expansion in \ as 

' 61 



X 



/y - 1 



(4.119) 



Using (|4.119jl in 1)4.118)1 and equating the various powers of y we get the following. 



The power of , y gives 



1 



1 ,D iC s 



(4.120) 



For?/ = — % we consider an expansion in \ as 



x 



+ a' + a[(y + i) + 



(4.121) 



Using (j4.121j) in ()4.118|) and equating the power of , \ 2 gives 



1 ,D zC\ 
1 ± -4/1 -4 

2 V v v' 



(4.122) 



To find the pole at infinity, we use the mapping V = \ Therefore equation ()4. 118)1 
transforms to 



x \ t} E)-t 2 X '(t,E) + - 



3 t 2 



1 t 2 



+ 



' [Et 2 — A — Bt — C f '' 



4(t 2 + l) 2 2t 2 + l 



v(f + 1 

For the point at infinity we expand \ as 



D- 



t 2 + 1 t 2 + 1 



(4.123) 



X (t, E) = d + dit + d 2 t 2 + 



Using ()4.124j) in 1)4.123)1 and equating the constant term we get 



(4.124) 



do = ±\/- 
V v 



(4.125) 
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and equating the term in t we have 

di = (4.126) 

Now equating the sum of residues due to fixed poles moving poles and that at 
infinity,to zero we have the following relation. 

b 1 + tf 1 +n = d 1 (4.127) 

Substituting the values of bi,b[ and d\ in the above relation yields the required 
condition given below 

B 4 + 16AV(n + l) 4 + 2AAB 2 u(n + l) 2 + AAB 2 D - AB 2 u(n + l) 2 

±32ABuVA^(n + l) 3 ± AABv^M>{n + 1) - AB 2 v - AA 2 u 2 (n + l) 2 
+ l6A 2 Du(n + l) 2 ± l6ABDVA^(n + 1) ± 8B 3 y/Au(n + 1) - AA 2 C 2 (4.128) 
This can be written in the compact form as: 

B±2(n + \)V^A\ A + A(AD - u) [b ± 2(n + l)V^Af - AA 2 C 2 = (4.129) 

The above condition agrees with those given in [9,10] 

4.8 V(x) = A cosh 2 ^/vx+B cosh y^x+C coth y / z/xcschy / z/x+ 
Dcsch yz/x 

The potential is given as 

V(x) = A cosh 2 \/~ux + £? cosh y^a; + C coth v^cschv^a? + -Dcsch 2 \/vx (4.130) 
We effect a change of variable by 

y = e^ x (4.131) 
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The QHJ equation is (h — 1 = 2m) 



p 2 (x, E) - ip'(x, E)-(E-V)=0 



(4.132) 



In the new variable the QHJ equation is 



Let 



p 2 (y,E)-iV^P'(y,E) 

E - A ^ 2 ^ 2 - B h {y2+lf - c i^ 



p(y, E ) = -iVvyfi 



-D 



% 2 



(y 2 - if 



(4.133) 



(4.134) 



In terms of the above transformation the QHJ becomes 



+ 



vy 



l 

"V 



l 



e - a-^ 2 + 1) 2 - B—( y 2 + 1) 2 - e rr* - £>- 



,2y(y 2 + l) _ 4y 2 



(</ 2 - l) 2 (y 2 - l)^ 



Let 



X = + 



2y 



= 
(4.135) 

(4.136) 



Therefore the above equation becomes 



.2 , / 



+ 



vy 



x +x+ ^ 



l 



.^V* +1) + -°(y 2 -i) 2 



2y(y 2 + 1) n 4y 



D 



(y 2 - if 



= 
(4.137) 



X has poles at y — and at y = ±1 and there are moving poles between the turning 
points. We assume that there are no more poles in the complex line. We find the 
residues for each of these pole below. 
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For y = we consider an expansion in x as 

X= -o H h cto + aiy H 

Using (j4.138|) in (|4.137jl and equating the power of A gives 



(4.138) 



bo = ±i 



(4.139) 



The power of p- gives 



6i = 1 + 



B 



For y = 1 we consider an expansion in x as 



+ a' + - 1) + 



Using (|4.141|) in (|4.137|) and equating the power of jzzff gives 



1 



,D 



For y = — 1 we consider an expansion in x as 



+ a + a/(j/ + 1) 



Using ()4.143|) in ()4.137|) and equating the power of , ^ 2 gives 



7 II 



1± \/l + 4( 



(4.140) 



(4.141) 



(4.142) 



(4.143) 



(4.144) 



To find the pole at infinity, we use the mapping y = \- Therefore equation ()4.137|) 
transforms to 
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we expand x(t,E) as 



X(t, E) =d + d x t + d 2 t 2 + ■ 



(4.146) 



Using ()4.14(jj) in ()4.145|) and equating the constant term we get 



M 

V Au 



Equating the term in t we have 



di 



B 
4udn 



(4.147) 



(4.148) 



Now equating the sum of residues due to fixed poles, the moving poles, and that at 
infinity, to zero we have the following relation. 



h + b[ + b'[ + 2n = d x 



(4.149) 



Substituting the values of b^b'^b'l and d\ in the above relation yields the required 
condition This can be written in the compact form as: 



B±2(n+l)VvA - A(4D + v) B±2(n + 1)ViM + AA 2 C 2 = (4.150) 



The above condition agrees with those given in [9,10] 



4.9 V(x) = Ae 2 ^ x + Be^ x + Ce~^ x + De~ 2 ^ x 



The potential is 



V(x) = Ae 2 ^ x + Be^ x + Ce~^ x + De~ 2 ^ x 



(4.151) 



The QHJ equation is given by (% — 1 = 2m) 



p 2 (x, E) - ip'(x, E)-[E- Ae 2 ^ x - Be^ x - Ce~^ x - De- 2 ^ x = (4.152) 
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We use a change of variable by 



y = e 



sfvx 



(4.153) 



Therefore the QHJ transforms to 



Let 



p 2 (y, E) - iV^yp'(y, E)-[E- Ay 2 - By - C- - d\] 

y y 



p(y, E ) = -iy/vy<i> 



= 



(4.154) 



(4.155) 



Therefore the above equation becomes 



1 1 



2y 



Ay 2 vy 1 



E — Ay 2 — By — C D— 

y y 



Let 



X = + 7T 



With this transformation we get 



2 /l 1 

X + X + H 9 



E-Ay 2 -By-C-- 

y y 



o 



(4.156) 



(4.157) 



(4.158) 



X has poles at y = and there are moving poles between the turning points. We 
assume that there are no more poles in the complex line other than a pole at infinity. 
We compute the residue for y = 

For y = we define 



b 2 bi 
X = -o H \-a + a 1 y + 

y y 



(4.159) 
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Substituting (|4.159j) in (|4.158j) we get 



H h a + 



+ 



) 2 + (-2^|-^ + a 1 + ...) + ^ 

y y % 



1 „ 1 



£ — Ay — By — C D 



y y 



Equating the coefficient of ^we get 



(4.160) 



bo = ±1 



(4.161) 



Equating the coefficient of p-we get 



&i = 1 



C 



(4.162) 



We assume that infinity is a pole and compute the residue at infinity, we which 
we use the mapping 



Therefore (|4.158|) transforms to 

X 2 (t,E)-t 2 x'(t,E) + ^ + - 

4 v 



E-A--B--Ct- Dt 2 
t 2 t 



(4.163) 



(4.164) 



We use a Laurent's expansion of the form 



X(t, E) = d + dit + d 2 t 2 + • • • (4.165) 
Using (|4.165|) in ()4.164|) and equating the of t the constant term yields 



do = ±\ 



and the coefficient of t yields 



dx 



B 
2d n v 



(4.166) 



(4.167) 
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Now equating the sum of residues due to fixed poles, the moving poles, and that 
at infinity, to zero we have the following relation. 

b 1 +n = d 1 (4.168) 

Substituting the values of b± and d\ in the above equation we get. 

±^L + (n + l) = ±^L (4.169) 

The above on simplification gives the desired condition given below which agrees 
with that given in [9,10]. 

2(n + l)y/vAD = ±BVT5 ± CVJ (4.170) 

4.10 Quartic Oscillator 

We end this chapter with a short analysis of quartic an-harmonic oscillator and give 
some remarks on polynomial potentials of degree different from six. 
First we consider the x 4 oscillator with 

V(x) = ax + (3x 2 + -fx 3 + 5x A , 5>0 (4.171) 

We ask whether this model is QES for any choice of parameters. We repeat the 
analysis given for sextic oscillator, assuming that the point at infinity is an isolated 
singular point, a pole of some order m. We therefore substitute 

p(x, E) = b m x m + b m ^x m ~ x + ■■■ (4.172) 

in the QHJ equation and determine the constants b m . For m > 2 we find that 
b m = and b 2 = ±iV6. Thus corresponding bound state wave-function for large x 
will behave as 

*p(x) ~ exp (ijx 3 ) (4.173) 
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and for both the choices ±iv^ for & 2 , one gets wave- function which grows either 
at + infinity or — infinity. Hence the assumptions, that the point at infinity is an 
isolated singular point of QMF, is inconsistent with QHJ for real parameter a, f3, 7 
and 5. Thus, x A oscillator does not lead to any choice of real parameter. However 
for complex parameters, one can get the known results for QES quartic model [18]. 

4.11 Summary and Observations 

From our study in this chapter we arrive at the following conclusions. 

1. For the QES models, the QMF corresponding to the algebraic part of the 
spectrum has singularity structure very similar to the exactly solvable models. 

2. The integer n appearing in the exact quantization condition is just the number 
of moving poles of QMF in the complex plane. In the case of exactly solv- 
able models the moving poles are in a one to one correspondence with the real 
nodes of the wave-functions, but a corresponding statement for QES model 
is not true. This result and some other interesting properties will be explic- 
itly demonstrated in the next chapter for bound state wave-functions of QES 
models. 

3. The integer n in the right hand side of the quantization condition appears in 
the condition of quasi-exact solvability of the potential. 

4. The condition of quasi-exact solvability is equivalent to our assumption about 
the behavior of QMF at infinity, reflecting the simplification of the singularity 
structure for the QES bound states. 



Chapter 5 



CALCULATION OF 
WAVE-FUNCTIONS FOR QES 
MODELS 

In the previous chapter we have seen that the condition for quasi-exact solvability 
arises from a simple requirement on the behavior of QMF at infinity. We continue our 
study of QES models and take up an investigation of the wave-functions. We find that 
the wave-functions can be computed by proceeding as in the case of exactly solvable 
models. We begin with our simplifying assumption mentioned in the previous chapter 
for the QES models and proceed in the same fashion as for the case of exactly 
solvable models in chapter 3. Thus the QMF is meromorphic and the corresponding 
residues at the poles are known, and also the behavior at infinity is known, with this 
information the bound state wave-functions can be obtained as in chapter 3. We 
give our results for two potential models viz. the sextic oscillator and the hyperbolic 
potential. This study reveals a new interesting feature of the zeros of the wave- 
functions, which will be discussed at the end of this chapter. 



63 



CHAPTER 5. CALCULATION OF WAVE-FUNCTIONS FOR QES MODELS 64 

5.1 Sextic Oscillator 

The potential for the sextic oscillator is: 

V(x) = ax 2 + px A + -fx 6 , 7>0 (5.1) 

with the following values for a, (3, 7 and the condition for fi,n,p where p stands for 
parity 

a = b 2 - a(3 + 2n), (3 = 2ab, 7 = a 2 , 4/i + 2p = 2n with p = or 1 
The QHJ equation is (7i = 1 = 2m) 

p(x, E) - ip'(x, E)-(E- ax 2 - (3x A - 7X 6 ) = (5.2) 

We assume that the point at infinity is a pole. Therefore p(x, E) behaves as x n for 
some n 

p(x, E) ~ x n 

for large x. Hence p(x, E) takes the form for large x. 

p(x, E) = CL3X 3 + a 2 x 2 + a\x + a + O(-) (5.3) 

x 

where do, a\ ■ • ■ , (13 are constants, on the assumption that p(x, E) have no other 
singular points and substitute (|5.3|) in equation (|5.2p . Next we equate the coefficient 
of powers of x 6 to zero, gives 

^ + 7 = (5.4) 

Since 7 = a 2 , we have 

03 = ±ia (5-5) 

As 03 has two values, the correct value is fixed by the condition of square inte- 
grability on the wave function. 
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if)(x) = exp \ i J p(x, E)dx S j = exp yi J (a 3 x 3 + aix)dxj 

If the above integral have to bounded at infinity, the we require that 

a 3 = +ia (5.6) 

Next equating the coefficient of successive powers x 5 , x 4 , . . . to zero we get 

a 2 = (5.7) 

a, = (5.8) 
a 3 

a = (5.9) 

Hence 

oi = ib (5.10) 

Therefore p(x,E) becomes 

n —i 

p(x, E) = h iax 3 + ibx + c (5-H) 

k=l X ~ Xk 

To determine x/-, or equivalently P(x) = Y\k=i( x ~ x k), we substitute (I5.11|) in ()5.2|1 

and get 

- [E - ax 2 - px 4 - 7X 6 ] = (5.12) 

Therefore, the above equation becomes 

pi pi pi pn 
c 2 + 2ax 3 — + Hbcx - 2i—c + 2bx— + 2iacx 3 \-b- E - 2anx 2 = (5.13) 

Equating the coefficient of x 3 term we have 

2iacx 3 = => c = 
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and hence we have 

pi pi pn 

2ax 3 — + 2bx— - — + b- E - 2anx 2 = (5.14) 

The above equation thus gives the following differential equation in P(x) 

P" - P'{2ax 3 + 2bx) — P(b — E — 2anx 2 ) = (5.15) 

We get the expression for energies and wave functions for various values of n as 
follows: 

We will derive explicit form of wave-functions for n = 0, 1 and 2. Later we will 
discuss the general form of the wave-function for arbitrary n. The general strategy 
for obtaining the wave-functions is the same as discussed for exactly solvable models 
in chapter 3. 

Wave-function for n=0: Only one energy level can be solved in this case. Since 
the number n, representing the number of moving poles is zero (j5.11j) . with c=0 as 
already found, becomes 

p(x, E) = tax + ibx (5.16) 
and hence the wave-function is given by 

ip(x) = exp J p(x)dx^j = exp J [iax 3 + ibx]dx^ = exp ^— a— — ^~2~^ 

and the corresponding energy is obtained from ()5.15|) by equating the constant term 
and is given as 

E = b. (5.18) 

Wave-function for n=l: In this case we take P(x,E) to be a first degree 
polynomial, (x — Xo). There ()5.15j) gives, xq = and the energy is given as 

E = 3b. (5.19) 
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Therefore the wave-function comes out to be 



ip(x) = Nx exp | — a— — b— j . 



(5.20) 



Wave-function for n=2: We seek a solution of (j5.15|) with P(x) as a second 
degree polynomial. Substituting P(x) as 

P(x) = ao + ol\x + a 2 x 2 . (5-21) 

Using the above equation in ()5.15j) and comparing different powers of x gives 

a x = (5.22) 

Aaa - a 2 (5b - E) = (5.23) 

a (b -E)- 2a 2 = 0. (5.24) 



The last two equation have non-trivial solution for ao and a 2 only if the 

4a (5b -E) 
'b-E) -2 

This gives two energy eigen- values 



(5.25) 



E = 3b ± 2Vb 2 + 2a. 



To get the wave function we compute ot\ and«2 from equation (|5.23J) and ()5.24|) and 



use 



/ \ . 2a 2 x . o 

= —i : + tax + ibx. 



ao + a 2 x 2 

Therefore the wave function is given by: 



(5.26) 



ip(x) = N exp (^i J p(x)dx^j = N exp (^i J [—i — -j— ^ — - + iax 3 + ibx\dx^j . (5.27) 



ip(x) = N(a + a 2 x ) exp I ~ b— 



(5.28) 
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where N is the normalizing factor. The value of «o is given by 

5b — E . 
«o = (5.29) 

Replacing the value of cto and energy value E in the above equation one gets the 
expression for wave-function as 



E — 56 / x 4 x 2 \ 
ij){x) = N b ± Vb 2 + lax 2 - 1 exp -a b— . (5.30) 



4a 



4 2 



The wave-functions and eigen-values explicitly obtained for the cases n = 0, 1 
and 2 agree with the known results [6]. 

For an arbitrary value of n the polynomial P(x) will be obtained by solving (j5.15|) . 
If we take P(x) to be of the form 

P(x) = J2a k x\ (5.31) 

k=0 

then the differential equation (|5.15|) leads to a set of homogenous equations for 
the corresponding coefficients ao, «i, • • • , ot n . These equations will have a non-trivial 
solution only if determinant of the coefficients vanishes. This condition will determine 
the energy eigen- value, corresponding to each eigen-value we can find the coefficients 
ao, a\, • ■ ■ , a n . Thus we get n independent wave-function each having the form 

ip(x) ~ P(x) exp f-a^- - b^j . (5.32) 

Notice that all these eigen-functions corresponding to a fixed value of n have a 
olynomial of the same degree n as a factor. Thus for a fixed value of n, and hence 
for a given set of potential parameters, wave-functions for all the states which can 
be solved have the same number of zeros equal to n. If these levels are arranged 
according to increasing energy, the number of zeros on the real axis (nodes) will 
increase. Hence the number of complex zero will decrease with increasing energy. 
This feature appears to be a general property of quasi-exactly solvable models. 
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5.2 Hyperbolic Potential 

The hyperbolic potential is 

A B 

V{x) = g- + o C cosh 2 x + D cosh 4 x, (5.33) 

cosh x sinh x 

where 

A = 4( ai -i)( ai -|), 
5 = 4( a2 -l)( S2 -|), 
C= [g 2 + 4g 1 (s 1 + s 2 + /i)], 
L> = g 2 . 

The QHJ equation is (ft = 1 = 2m) 

A B 

p 2 (x, E) — ip'(x, E) — [E + — -g -5- + C cosh 2 x - D cosh 4 x] =0. (5.34) 

cosh x sinh a; 

We effect a transformation by 

y = cosh a; (5.35) 

The QHJ equation in the new variable is 

A B 



p(y, E) - iTi^y* - lp'(y, E)-[E + - 2 ~ + Cy 2 - Dy 4 } = (5.36) 

V yZ yi _ j_ 

Let 

p(y, E) = -iy/y* - l^y, E) (5.37) 
Therefore the QHJ equation becomes 



1 y 



2y 2 -l 
Let 



I(^ + ^ + ^^ + 7-^T + ^ 2 -^-°- ^ 



i v 

2y 2 -l 



X = 0+o7I^T- ( 5 - 39 ) 
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Therefore the above equation becomes 

X 2 + x' + l T ^ V2 -l^ + ^[E + - 2 -^- + Cy"-Dy'] = (5.40) 
A[y z — l) z 2y z — l y —1 y y — I 

X has poles at y — 0, and y — ±1, and there are moving poles between the turning 

points. We assume that there are no more poles in the complex plane other than 

a pole at infinity. We will first compute the residues at y — 0, ±1 and then in the 

general form of \ (I5.51j) the constants b\, b\ and b'[ will be known and then we give 

the general form of the wave-function. 

Computation of residues: For y = let 

x = — + a + a 1 y^ (5.41) 

y 

Therefore equation (J5.40j) becomes 

r b\ l9 r bi . 3 y 2 11 

[- + a + aiy + ■ ■ f + [ — r + ai + •••] + -- 



y ' ^ ' — ' J ' L y* ■ - ■ J ■ 4(^-1)2 2 y«-l 

+ + ^ - -At + ^ 2 - = ( 5 - 42 ) 

Equating the coefficient of on both sides gives 

&i = ^[l±(4 Sl -2)] (5.43) 

By the condition of square integrability, the positive sign has to be taken. Hence the 
value of b\ is 

h = 2 Sl - \ (5.44) 

For y — 1 let 

b' 

X= -^— + a' + a' 1 {y-l) + --- (5.45) 

y- 1 

Therefore equation ()5.40j) becomes 

[-^- + a' + a' 1 (y-l) + ...] 2 + [--i- + a' 1 + ...] + 3 ^ 1 1 



(y-l)a 1 J 4(y 2 -l) 2 2 y2 
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+^-^[E + - 2 -^— + Cy 2 -Dy 4 ]=0 

y — i y y i 

Equating the coefficient of ( y - 1 1 ) 2 on both sides gives 

6' 1 = 1[1±(2 S2 -1)] (5.46) 

By the condition of square integrability, the positive sign has to be taken. Hence the 
value of b[ is 

b[ = s 2 (5.47) 



X = + < + <(y + !) + •■■ (5.48) 



For y — — 1 let 



Therefore equation ()5.40|) becomes 
[^ T + < + <( 2 / + l) + ---] 2 + h— ^ + < + "--] + 3 //J 1 1 



+ l ' u ' 1V * ' ' ' J L (y + l) 2 J 4(y 2 -l) 2 2y 2 -l 

+^ i T^ + 4-^T + ^ 2 -^ 4 ] = ° 
y — l y y — l 

Equating the coefficient of ( - y + 1 1 ) 2 on both sides gives 

U{ = ± (2s 2 - 1) ] (5.49) 

By the condition of square integrability, the positive sign has to be taken. Hence the 
value of b" is 

b'[ = 8 2 . (5.50) 

For the fixed poles at y — 0, ±1 let % have the form 

P\y) h b[ b'l . . 

x = -Err + - + — ~t + ~~t + ay + 02, 5.51 
P{y) y y - 1 y + 1 

where ci and c 2 are constants to be determined. This form of x is because the 
equation ()5.4()|) has a y 2 term. 
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Form of wave- function: Using (|5.51|) . equation (|5.40|) transforms to 



p'{v) , 61 , b 2 



n 2 



P(y) y y - i y + 1 
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+ c x y + c 2 



+ 



^(y) f^(y)V 
^(y) V ^(y) / 



y 2 (y _ 1)2 (y + l) 2 



+ Ci + - 



1 1 



1 m -4 

y - 1 r 



y 2 -i 

and for large y equating the coefficients of y 2 to zero gives 

c\ - D = 0, 



4(y 2 -l) 2 2y 2 -l 
Q/ 2 - Dy 4 } = 



(5.52) 



hence 



c x = iv^D = ± 9l . 



For large y equating the coefficients of y gives 



2c lC2 = 0, 



(5.53) 



hence 



c 2 = 0. 



(5.54) 



The correct sign of c\ is fixed by square integrability and is found to be c\ = —q\. 
As the potential is symmetric, there are moving poles on either side and hence we 
take P(y) to have the form given below 



P(y) = Il(y 2 -yl)- 

fc=i 

The wave-function for this model is computed as follows 

if)(y) = exp [ [x- \ 2 V ]dy. 



(5.55) 



ly 1 



exp 



rP'{y) , h b[ 



P{y) y y - 1 y + 1 



(5.56) 
(5.57) 
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On integrating and substituting the values of bi, b 2 and c\ we get the expression for 
the wave-function in the x variable as 

tp(x) = (cosh 2 (sinh 2 x) S2 ~ 3 exp ( — cosh 2 x ) (cosh 2 x — y 2 ). (5.58) 

^ 2 ' k=i 

Computation of energy-eigenvalue: We shall now show how our analysis 
leads to the correct answer for energy spectrum. With c 2 = and substituting the 
values of b\, b 2 a.iad c\ (|5.52j) takes the form 



P"{y) P'{v) 
P(y) P(y) 



22 4SiS 2 4sis 2 s 2 

+mv + — - — + 



(45i-l) | 2s 2 



y 

S-2 



2s 2 

y-1 ' y + 1 



2qiy 



+ 



s 2 



y-1 y+1 y+1 y-1 y-1 y+1 



So y y 

2s 2 qx — — -2s 2 <?i 



- 4«igi - -- 



1 1 



E + A B 1 



2y 2 



c- 



D- 



y 



y+1 

l 



y 2 -v 



y-i 

(5.59) 



Using (|5.55|) in ()5.59j) we get the following equation. 

v 2y 2 

2-^ „,2 _ „,2 



E 

fc=i 



i i 

+ 



(y + yk) 2 (y - ykY 



+ E 



, fc =i y - y% 
2y 



*=i y 2 -ylv y 



(4 Sl - i; 



2s, 



2s, 



2 2, 4s l S 2 4siS 2 . S 2 

+[qfy + — — + 



S-2 



+ 



y-1 y+1 

y 



2q\y 



y-1 y+1 y+1 y-1 y-1 y+1 



2s 2 q x 



E + A B 1 
+ 



y 2 -1 y 2 - 1 2 y 2 — 1 



C- 



D- 



y 



y + 1 

A 



2s 2 qi 



y 



y-i 

(5.60) 



Multiplying the above throughout by i and integrating along a closed contour en- 
closing y = we get the expression for energy as 



E = -4 



n 2 



si + s - 2 



8si 



qi 



+ E 

k=l 



y k z 



(5.61) 



Using y\ = ^ the above equations for energy become 
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E = -4 



i 2 



si + s - 2 - 



Changing to = £ and j/j^ = (|5.6U|) becomes 

n. -i "I 2 n 



E 



2E 

fc=l L 

1 



£-1 



(e - ao 2 

2s 2 ' 



+ 2 



+ H 



(4 Sl - 1) + 4s 2 
1 



. 2 fc=i ?* 



(5.62) 



- 4s 2 gi 



£-1 



+ ^ r^ + — T - + C- 



D 



e 2 



- 2gi£ 

4s x gi - 
= 



E 

fc=l L 



i i 



(5.63) 

e-i 2 e- 1 e-i e-i 

and integrating ()5.63|) over £ around a closed contour, enclosing only one of the 
points £j and repeating for (i — 1, 2, • • ■ , n) we get the following result 



n i 

n 



£i s 2 



k=l S» sfc si si 1 2 

i = 1, 2, • • • n 



(5.64) 



The results fl5.58j) . ()5.62|1 and ()5.64j) are in agreement with those given in [6] 

The general feature of the zeros of the wave-function for the sextic oscillator are 
also true for the QES hyperbolic potential. In particular, for a given potential it 
is correct that all the exactly solvable wave-functions have the same total number, 
(real and complex) of zeros. This feature is found to be correct for all QES model 
studied so far including the QES periodic potentials [12]. 

5.3 Concluding Remarks 

Our study of bound state wave-functions in this chapter shows the following similar- 
ities and differences between the exactly solvable and QES models. 

1. In both the models, the "QMF" turns out to be a rational function after a 
suitable change of variables. 
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2. In both the cases, the integer n in the right hand side of quantization condition 
coincides with the number of moving poles. 

3. For every bound state in one dimension the k th excited state wave-function 
have A;-nodes on the real axis. This statement is a general one and is true for 
all models including exactly solvable and QES potentials. The study in chapter 
3 shows that QMF for exactly-solvable models has moving poles which are in 
correspondence with the nodes of the wave-function. There are no poles off 
the real axis. However this property fails to be true for QES potentials where 
the QMF has poles off the real axis, in addition to the poles on the real axis 
corresponding to the nodes of the wave-function. 

4. For the QES potentials only a part of the energy spectrum and the correspond- 
ing wave-functions can be computed exactly. An interesting property of the 
QMF for all these levels is that the total number of moving poles is the same 
and equal to the integer n of the quantization. 

5. Different values of integer n correspond to different QES potentials within a 
family, and it does not refer to different excited state of a single potential, as 
was the case for exactly-solvable model. 



Chapter 6 

CONCLUSIONS AND 
OUTLOOK 

In this thesis, we have studied exactly solvable and QES potentials in one dimensional 
quantum mechanics with in the frame work of QHJ formalism. The following results 
have been obtained. 

1. The eigen- values and eigen-functions of the exactly solvable models can be 
obtained by a very simple and elegant method, which makes use of elementary 
results from theory of complex variables. 

2. The non-trivial input in this analysis is the singularity structure of the QMF. 
Besides this, a change of variable is needed to transform the QHJ equation to 
a Riccati form with rational functions as coefficients. Having done this, it is 
very easy to identify the fixed poles and the corresponding residues. It is the 
location and the number of moving poles which present some difficulty. For 
a large class of exactly solvable and QES models studied by us, the number 
of moving poles for the bound states turns out to be finite. In addition, the 
behavior of QMF for large values of independent variables has been very simple 
to read from QHJ equation. All these observations can be summarized in one 
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sentence by saying that the QMF is a rational function after a suitable change 
of variables. 

3. The quantization condition as given by Leacock and Padgett, is applicable to 
separate systems which can be reduced to one dimensional problems. It will 
be interesting to formulate an exact quantization condition for non-separable 
systems in higher dimensions, and investigate its relation to the well known 
existing semi-classical schemes and to see applications to chaotic systems. 

4. For other models which are not exactly solvable, one has to device and ap- 
proximation scheme. Here again some idea about the knowledge of location of 
moving poles has any relation to classical trajectories. 

5. We have tried to study the QHJ formalism for an-harmonic oscillator which 
is a test case of any computational scheme. One can compute the asymptotic 
value of the QMF for large x and one can use this answer as an input for nu- 
merical integration of QHJ equation. Detailed investigation is in progress and 
interesting approximation scheme for an-harmonic oscillator is expected from 
the preliminary results. The simple form of QHJ equation offers a possibility 
of several analytic approximation schemes also. 

6. When computing bound states for some potentials such as Rosen Morse hy- 
perbolic potential, it is found that applying boundary condition p(x, E) ^5 
p c i(x,E), carefully leads one to select different residues for different ranges of 
potential parameters. Excepting this and proceeding, further analysis leads to 
different set of energy spectrum and wave-functions for such different ranges of 
parameters in the potential. This is consistent with the known result on phases 
of super-symmetry in Rosen Morse potential [13]. Other such potentials, for 
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example trigonometric Scarf potential [14], which exhibit different phases for 
different ranges of potential parameters, can also be investigated within our 
frame work. 

7. It must be remarked that, the QHJ formalism as presented in this thesis is 
applicable to bound states only. Modifications will be needed to apply this 
formalism to continuous energy solutions. The requirements such as the quan- 
tization rule, square integrability etc., are no longer applicable. For such cases 
we must accept all possible combination of residues, consistent with the other 
equations of the theory, and proceed to analyse the consequences. In fact 
analysis of this type has been performed for some of the QES and exactly solv- 
able periodic potentials, and QHJ formalism leads to the full set of band-edge 
wave-function and corresponding energy eigen- values [12,15]. 

8. The QHJ formalism offers advantages from the pedagogical point of view. The 
understanding of method and results requires only the basic understanding of 
the theory of complex variables. In this connection we mention periodic poten- 
tials, both exactly solvable and QES, and PT symmetric complex potentials 
[16] which can be handled with equal ease within the QHJ approach [12,15, 
17]. 
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